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ABSTRACT 

This  paper  demonstrates  the  effectiveness  of  the  zero-forcing,  minimum  mean-squared  error,  and 
constant-modulus  equalizers  in  improving  the  BER  performance  of  iNET-formatted  SOQPSK-TG. 
The  equalization  algorithms  leverage  the  existence  of  known  bit  sequences  in  the  preamble  and 
ASM  fields  of  the  iNET  packet  to  realize  data-aided  equalizers.  The  effectiveness  of  these  equal¬ 
ization  techniques  over  ten  test  channels,  derived  from  channel  sounding  experiments  at  Edwards 
AFB,  was  evaluated.  The  BER  curves  for  nine  out  of  the  ten  test  channels  display  the  desirable 
“waterfall”  shape.  The  BER  curve  for  the  remaining  test  channel  displays  a  “BER  floor.”  Fortu¬ 
nately,  the  “BER  floor”  is  low  enough  to  allow  error  correcting  codes,  such  as  the  iNET  LDPC 
code,  to  correct  the  errors  and  provide  virtually  error-free  performance. 


INTRODUCTION 

Multipath  interference  occurs  when  the  radio  signal  form  the  airborne  transmitted  arrives  at  a 
ground-based  receiver  via  multiple  propagation  paths.  Usually,  one  of  the  propagation  paths  is 
the  line-of-sight  propagation  path  whereas  the  others  are  due  to  reflections.  Multipath  interference 
continues  to  be  the  dominant  cause  of  link  outages  in  aeronautical  telemetry. 

Equalization,  using  the  constant  modulus  algorithm  (CMA),  in  the  context  of  aeronautical  teleme¬ 
try  has  been  investigated  before  [1,  2,  3,  4].  However,  laboratory  tests  with  hardware  prototypes 
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Figure  1:  The  iNET  packet  structure  used  in  this  paper. 


produced  less  than  compelling  results  [4];  these  results  have  raised  questions  on  the  applicability 
to  aeronautical  telemetry  of  purely  blind  equalization  algorithms. 

Here  we  investigate  a  data-aided  approach  to  equalization.  In  data-aided  equalization,  the  equalizer 
filter  coefficients  may  be  computed  from  the  multipath  channel  coefficients.  The  question  is,  “how 
does  the  receiver  know  what  the  channel  impulse  response  is?” 

The  answer  to  this  question  lies  in  the  iNET  packet  structure.  iNET-formatted  transmissions  in¬ 
clude  a  128-bit  preamble  and  64-bit  attached  sync  marker  (ASM)  preceding  a  block  of  data  bits 
(at  least  6144  bits:  an  LDPC  codeword):  see  Figure  1.  Because  the  preamble  and  ASM  bits  are 
known,  the  receiver  can  compare  the  received  signal  to  a  locally  stored  copy  of  the  SOQPSK-TG 
signal  corresponding  to  the  preamble  and  ASM  bit  fields.  This  comparison  is  capable  of  producing 
estimates  of  the  frequency  offset,  noise  variance,  and  multipath  channel  coefficients  [5].  And  from 
the  multipath  channel  coefficient  estimates,  the  equalizer  filter  coefficients  may  be  obtained. 

In  this  paper  we  evaluate  the  performance  of  three  popular  equalization  algorithms:  the  zero¬ 
forcing  (ZF)  equalizer,  the  minimum  mean-squared  error  (MMSE)  equalizer,  and  the  CMA  equal¬ 
izer.  The  ZF  and  MMSE  equalizers  are  data-aided  equalizers.  Traditionally,  the  CMA  equalizer  is 
a  completely  blind  equalizer;  but  here  we  use  the  MMSE  equalizer  filter  coefficients  to  initialize 
the  CMA  equalizer.  Relative  to  the  traditional  “center-tap  initialization,”  the  MMSE-initialized 
CMA  equalizer  as  a  shorter  convergence  time  and  lower  bit  error  rate. 


SYSTEM-LEVEL  DESCRIPTION 

The  bit  sequence  for  iNET  is  organized  as  outlined  in  Figure  1.  The  preamble  sequence  (PRE) 
is  CD98hex  repeated  eight  times  [6,  p.  48].  The  preamble  field  is  followed  by  the  attached  sync 
marker  (ASM)  field  defined  as  034776C7272895B0hex.  The  DATA  field  is  6144  randomized  data 
bits.  (These  bits  correspond  to  a  single  LDPC  codeword  in  the  coded  system.  Here,  we  evaluate 
the  uncoded  bit  error  rate  after  equalization.) 

The  transmitted  signal  is  SOQPSK-TG  whose  input  bit  stream  is  summarized  in  Figure  1.  The 
signal  propagates  through  a  frequency  selective  channel  and  experiences  a  frequency  offset  as  well 
as  the  addition  of  additive  white  Gaussian  noise. 

The  received  signal  is  filtered,  down-converted  to  I/Q  baseband,  and  sampled  (not  necessarily  in 
that  order)  using  standard  techniques.  The  focus  of  this  paper  is  on  the  application  of  equalization 
techniques  to  the  EQ  baseband  samples.  Because  SOQPSK-TG  is  a  nonlinear  modulation,  the 
equalizer  cannot  operate  on  the  symbols  in  the  same  way  it  does  for  linear  modulation  (cf.,  [7, 
Chapter  9]).  Consequently,  the  equalizer  must  operate  on  the  samples  of  SOQPSK-TG,  similar  to 
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Figure  2:  The  data  packet  format  and  high-level  signal  processing  explored  in  this  paper. 


the  way  fractionally  spaced  equalizers  operate. 

Because  the  preamble  and  ASM  bits  are  known,  the  samples  corresponding  to  the  preamble  and 
ASM  bits  are  used  to  estimate  the  frequency  offset,  channel  impulse  response,  and,  for  the  MMSE 
equalizer,  the  noise  variance.  Before  these  tasks  can  be  accomplished,  the  start  of  the  samples 
corresponding  to  the  preamble  bits  in  the  received  signal  must  be  known.  This  is  accomplished 
by  the  preamble  detector  block,  whose  algorithm  is  based  on  the  detection  algorithm  described  in 
[8].  Once  the  start  of  the  preamble  is  known,  the  frequency  offset  is  estimated  using  the  algorithms 
described  in  [5].  The  frequency  offset  is  used  with  a  complex-exponential  to  derorate  the  received 
data  to  remove  the  frequency  offset.  The  derorated  data  rd{n)  are  used  to  estimate  the  channel 
and  noise  variance  as  described  in  [5].  The  channel  estimates  h(n),  for  — Ar,  <  n  <  N2,  are  used 
to  compute  the  ZF  and  MMSE  equalizer  filter  coefficients.  For  the  CMA  equalizer,  the  channel 
estimates  h(n)  are  used  to  initialize  the  adaptive  filter  whose  update  is  based  on  the  CMA. 


THE  EQUALIZATION  ALGORITHMS 

The  equalizers  operate  in  the  system  configuration  shown  in  Figure  3.  [cf.,  Figure  2].  Here,  the 
derotated  samples  rain)  are  equalized  using  a  length  L\  +  L2  +  1  FIR  filter  defined  by  the  impulse 
response  c(n)  for  —Li  <  n  <  L2  to  produce  the  output 

^2 

y(n)  —  c(m)rd(n  —  m).  (1) 

m=—Li 

The  equalizer  output  forms  the  input  to  the  well-known  symbol-by-symbol  SOQPSK  detector  com¬ 
prising  a  detection  filter  operating  at  N  =  Tb/T  samples/bit  and  a  decision  process,  operating  on 
the  decision  variable  u(k )  at  1  sample/bit.  This  detector,  based  on  an  offset  QPSK  approximation 
of  SOQPSK-TG,  is  described  in  more  detail  in  [9,  10].  The  detectors  of  Figure  3  also  include 
a  phase  lock  loop  (PFF).  The  PFF  is  required  to  track  out  any  residual  phase  increments  due  to 
frequency  offset  estimation  errors.  A  timing  loop  is  not  required  because  timing  offsets  are  part  of 
the  channel  estimate  h(n).  In  what  follows,  we  will  organize  the  equalizer  filter  coefficients  into 
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derotated  data 
samples 
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diagram  for  the  CMA  equalizer. 


(L i  +  L2  +  1)  x  1  vectors  as  follows: 


Czf(— Li) 

Cmmse(— Li) 

'4ma(-£0' 

CZF  ~ 

Czf(O) 

,  CmmSE  — 

cmmse(0) 

( p ) 

1  CMMSE  ~ 

ccma(0) 

.  Czf(-^2)  . 

.  Cmmse (L2)  _ 

.  CCMa(L2)  . 

The  superscript  “(p)”  for  the  CMA  equalizer  coefficients  is  the  update  index  and  is  explained 
below. 

The  ZF  equalizer  is  a  filter  that  is  the  best  length-  (Li  +  L2  +  1)  FIR  approximation  to  the  “inverse” 
of  the  channel.  That  is,  c(n)  is  chosen  so  that 


h(n)  *  c(n)  ~  S(n  —  no)  (3) 

for  —L\  <  n  <  L2  where  “best”  minimizes  the  least  squares  error.  The  vector  of  filter  coefficients 
defined  by  this  criterion  is  given  by  [11] 


cZF=(HtH)  1  H^uno 


(4) 


where 


u 


no 


"O' 

0 

1 

0 

0 


^  n0  —  1  zeros 

|  Ni  +  N2  +  Li  +  L2  -  n0  +  1  zeros 


and  H  is  the  (N{  +  N2  +  L±  +  L2  +  1)  x  (L i  +  L2  +  1)  convolution  matrix  formed  by  the  estimates 
of  the  channel  impulse  response: 


hi-Nj 

hi-Ni  +  1)  hi-Ni) 

h(N2 )  h(N2  -  1) 

h(-Ni) 

h(N2) 

hi-Nt  +  1) 

h(N2) 

The  minimum  mean-squared  error  (MMSE)  equalizer  is  a  filter  that  minimizes  the  mean  squared 
error 


£ 


E 


|  s(n)  —  r(n )  *  c(n) 


(5) 
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The  solution  is  given  by  a  form  of  the  Wiener-Hopf  equations  [11].  The  (L i  +  L2  +  1)  x  1  vector 
of  MMSE  equalizer  filter  coefficients  are 


cmmse  — 


GGt  +  °^I 


3  2*-Li+l2+i 


G  is  the  T  L2  T  1)  x  (Ay  T  N2  T  Ti  T  L2  T  1)  matrix 


h(N2)  •••  hi-Nj 

h(N2)  •••  hi-Nt) 

h(N2)  •••  h(-Ni), 


(6) 


g  is  the  1  x  ( Li  +  L2  +  1)  vector  given  by 

g  =[h{Li)  K-L2)\ 


where  it  is  understood  that  h(n )  =  0  for  n  <  —Ay  or  n  >  Ay  (how  many  zeros  need  to  be 
prepended  and  appended  depends  on  the  relationship  between  L1  and  Ay  and  the  relationship  be¬ 
tween  L2  and  A' 1 );  and  a2s  is  the  signal  power.  This  solution  is  based  on  the  assumption  that 
SOQPSK-TG  samples  corresponding  to  a  sample  rate  of  2  samples/bit  are  approximately  uncorre¬ 
lated.  The  impact  of  this  approximation  is  discussed  in  [12]. 

The  CMA  equalizer  is  outlined  in  Figure  3  (b).  As  with  the  ZF  and  MMSE  equalizers  of  Figure  3 
(a),  the  lcngth-fL]  +  L2  +  1)  FIR  equalizer  filter  precedes  the  SOQPSK-TG  symbol-by-symbol 
detector.  Unlike  the  ZF  and  MMSE  equalizers  of  Figure  3  (a),  the  CMA  equalizer  is  an  adaptive 
filter  whose  coefficients  minimize  the  cost  function  Jcma  given  by 


Jc ma  (y(n)^j  =  E|(||/(n)|2  -i?2)2} 


(7) 


where 


(8) 


The  CMA  cost  function  measures  the  departure  of  the  equalizer  output  from  a  circle  of  radius  R2, 
but  does  so  without  any  phase  information. 


Using  the  steepest  descent  algorithm  to  drive  the  adaptation,  which  occurs  once  per  packet  here, 
the  ( L 1  +  L  )  T  1)  x  1  vector  of  filter  coefficients  for  packet  p  +  1  is  an  updated  version  of  the  filter 
coefficients  used  for  packet  p: 

ccma^  =  ccma  —  /^V  Jcma  (9) 

where  p  >  0  is  the  step  size  and  where 


V  Jcma  —  E  2 


y(n)y* 


n) 


R-l 


y{n ) 


r  (n, 


(10) 
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Figure  4:  A  block  diagram  of  the  simulation  procedure. 


with 


r  n  = 


rd(n  +  Li) 


r*dH 


(ii) 


L r*d(n  ~  L2) J 

A  critical  issue  for  any  adaptive  filter  is  convergence:  both  how  long  and  to  what  state.  Here,  we 
initialize  the  CMA  equalizer  filter  with  the  MMSE  filter  coefficients.  That  is 


c 


(0)  _ 

CMA  ~  CMMSE 


(12) 


given  by  (5).  Not  only  does  this  reduce  convergence  time,  but  also  it  helps  the  adaptive  filter  con¬ 
verge  to  a  state  the  produces  a  lower  bit  error  rate  than  that  achievable  with  center-tap-initialization. 


PERFORMANCE  RESULTS 

The  bit  error  rate  (BER)  performance  of  the  three  equalization  techniques  was  assessed  using  the 
simulation  environment  outlined  in  Figure  4.  The  simulation  parameters  were  the  following: 

1.  The  payload  data  rate  was  equivalent  to  10  Mbits/s  (the  equivalent  “over-the-air”  bit  rate  was 
10.3125  Mbits/s).  The  iNET-formatted  SOQPSK-TG  signal  and  channel  were  generated  at 
an  equivalent  sample  rate  of  2  samples/bit. 

2.  Because  the  channel  estimator  does  not  know  the  true  length  of  the  channel,  the  estimator 
used  values  for  Ni  and  N2  larger  than  any  of  the  test  channels.  These  values  were  Af  =  12 
and  N2  =  25  samples. 

3.  The  equalizers  used  L\  =  4  x  Ni  =  48  samples  and  L2  =  4  x  N2  =  100  samples.  Thus  the 
length  of  equalizer  filter  was  149  samples. 

4.  The  CMA  equalizer  filter  is  an  adaptive  filter.  After  initializing  the  filter  with  the  MMSE 
equalizer  filter  coefficients,  the  CMA  filter  was  allowed  to  converge  for  100  packets  before 
errors  were  counted. 
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5 .  The  simulations  were  performed  over  10  representative  channels  derived  from  channel  sound¬ 
ing  measurements  conducted  at  Edwards  AFB  under  the  M4A  program  [13].  The  test  chan¬ 
nels  are  summarized  in  Table  1  and  the  corresponding  frequency-domain  plots  are  shown  in 
Figure  5. 

The  simulated  BER  performance  is  shown  in  Figures  6  -  10.  In  all  cases  we  observe  that  the  ZF 
equalizer  has  the  worst  performance  and  that  the  MMSE  and  CMA  equalizer  have  about  the  same 
performance.  (The  lone  exception  is  the  simulation  results  for  test  channel  9  shown  in  Figure  10 
where  for  high  Eb/N0,  all  three  equalizer  achieve  the  same  performance.  This  is  to  be  expected 
because  test  channel  9  is  rather  benign  [see  Figure  5].)  The  fact  that  the  ZF  equalizer  has  the  worst 
performance  is  not  surprising.  The  ZF  equalizer  simply  “inverts  the  channel”  (this  phrase  derives 
from  the  frequency  domain  point  of  view).  For  channels  with  nulls,  the  “inversion”  restores  the 
frequency  content  of  the  desired  signal  in  the  frequency  band  surrounding  the  null.  This  restoration 
also  amplifies  the  noise  power  in  the  same  frequency  band.  The  end  result  is  a  phenomenon 
known  as  “noise  amplification:”  the  distortion  due  to  the  multipath  channel  is  corrected,  but  the 
signal-to-noise  ratio  (greatly)  decreased.  In  contrast  the  MMSE  and  CMA  equalizers  take  a  more 
measured  approach  to  “channel  inversion”  and  balance  the  impact  of  residual  multipath  distortion 
and  amplified  noise  on  their  respective  cost  functions  [mean  squared  error  (4)  for  the  MMSE 
equalizer  and  JCMA  (6)  for  the  CMA  equalizer]. 

Another  interesting  feature  of  the  simulation  results  is  that  all  test  channels,  except  test  channel  2, 
display  the  “waterfall”  shape;  a  large  decrease  in  BER  as  Eb/N0  increases  above  a  certain  value. 
Systems  characterized  by  BER  vs.  Eb/N0  curves  with  the  “waterfall”  shape  are  those  for  which 
equalization  and  additional  link  margin  can  achieve  a  desired  reliability  level.  The  good  news 
here  is  that  for  nine  out  of  the  ten  cases,  reliable  communications  is  possible.  Test  channel  2  is  a 
counter  example.  Here,  the  BER  performance  seems  to  “flatten  out”  as  Eb/N0  increases.  That  is, 
the  BER  curve  is  characterized  by  a  “BER  floor”  rather  than  the  “waterfall”  shape.  A  “BER  floor” 
is  undesirable  because  it  means  that  no  matter  how  much  Eb/N0  is  increased,  the  system  designer 
cannot  drive  the  BER  below  the  floor.  All  is  not  lost,  however.  If  the  “BER  floor”  is  low  enough 
(as  it  is  in  this  case),  then  a  good  error  correcting  code  (such  as  the  LDPC  code  defined  in  the  iNET 
standard)  can  remove  the  errors  to  produce  a  close-to-error-free  link. 


CONCLUSIONS 

This  paper  has  demonstrated  the  effectiveness  of  the  ZF,  MMSE,  and  CMA  equalizers  with  iNET- 
formatted  SOQPSK-TG.  The  equalization  algorithms  leverage  the  existence  of  known  bit  se¬ 
quences  in  the  preamble  and  ASM  fields  of  the  iNET  packet  to  realize  data-aided  equalizers. 
The  effectiveness  of  these  equalization  techniques  over  ten  test  channels,  derived  from  channel 
sounding  experiments  at  Edwards  AFB,  was  evaluated.  The  BER  curves  for  test  channels  1,  3-10 
display  the  desirable  “waterfall”  shape,  whereas  the  BER  curve  for  test  channel  2  displays  a  “BER 
floor.”  Fortunately,  the  “BER  floor”  is  low  enough  to  allow  error  correcting  codes,  such  as  the 
iNET  LDPC  code  to  correct  the  errors  and  provide  virtually  error-free  performance. 


8 


Figure  5:  Frequency-domain  plots  of  the  example  channels  from  channel  sounding  experiments  at 
Edwards  AFB.  In  each  plot,  the  thick  line  is  the  channel  frequency  response  and  the  thin  line  is  the 
power  spectral  density  of  SOQPSK-TG  operating  at  10.3125  Mbits/s. 
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Table  1:  Description  of  the  ten  test  channels  used  in  the  simulations. 


channel 

N, 

n2 

length 

environment 

01 

2  1 

7 

9 

Taxiway  E 

02 

2 

17 

20 

Taxiway  E 

03 

1 

22 

24 

Taxiway  E 

04 

5 

13 

19 

Takeoff  on  22L 

05 

1 

1 

3 

Cords  Road 

06 

1 

2 

4 

Cords  Road 

07 

0 

4 

5 

Cords  Road 

08 

2 

3 

6 

Black  Mountain 

09 

1 

1 

3 

Black  Mountain 

10 

2 

3 

6 

Land  on  22L 

Figure  6:  Simulation  results  for  test  channel  1  (left)  and  test  channel  2  (right). 
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Eb/NO  (dB)  Eb/NO  (dB) 


Figure  7:  Simulation  results  for  test  channel  3  (left)  and  test  channel  4  (right). 
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Figure  8:  Simulation  results  for  test  channel  5  (left)  and  test  channel  6  (right). 
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Figure  9:  Simulation  results  for  test  channel  7  (left)  and  test  channel  8  (right). 


Figure  10:  Simulation  results  for  test  channel  9  (left)  and  test  channel  10  (right). 
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